arXiv:1507.03595vl [cond-mat.mes-hall] 13 Jul 2015 


Channel Blockade in a Two-Path Triple-Quantum-Dot System 

M. Kotzian/’[^ F. Gallego-Marcos,^’^ G. Platero,^ and R. J. Haug^ 

^ Institut fiir Festkorperphysik, Leibniz Universitdt Hannover, Appelstrasse 2, 30167 Hannover, Germany 
^Instituto de Cieneia de Materiales, CSIC, Cantoblaneo, 28049 Madrid, Spain 

(Dated: July 15, 2015) 

Electronic transport through a two-path triple-quantum-dot system with two source leads and one 
drain is studied. By separating the conductance of the two double dot paths, we are able to observe 
double dot and triple dot physics in transport and study the interaction between the paths. We 
observe channel blockade as a result of inter-channel Coulomb interaction. The experimental results 
are understood with the help of a theoretical model which calculates the parameters of the system, 
the stability regions of each state and the full dynamical transport in the triple dot resonances. 

PACS numbers: 73.21.La, 73.23.Hk, 73.63.Kv, 85.35.Ds, 85.35.Gv 


Triple quantum dots (TQDs), which have been imple¬ 
mented only recently offer the possibility of ana¬ 

lyzing new fascinating properties which are not present 
in double-quantum-dot systems. These new properties, 
to name a few, include interference phenomena between 
different transport channels giving rise to dark states in 
triangular [SHE] and linear [9] dot distributions and long 
distant coherent states in TQDs [9UT3]. TQDs are, as the 
smallest qubit chain, a step towards more complex archi¬ 
tectures needed in quantum computation. They allow 
for novel applications in the field of quantum informa¬ 
tion processing, like for example as exchange-controlled 
spin qubits [TH [15] or as current rectifiers [TJ |T6] . They 
provide as well the implementation of quantum cellular 
automata processes, a combination of charging and re¬ 
configuration events in the system being a crucial pro¬ 
cess in quantum information [niiiH]. Coherent electron 
transfer using adiabatic passage was proposed for TQDs 
in series [T9|. Furthermore, decoherence due to charge 
fluctuations is reduced in a TQD-based coded qubit as it 
involves a decoherence free subspace [I5l|20]. 

Our system is a triangular-shaped TQD with one lead 
attached to each dot thus consisting of two double-dot 
paths. A triangular geometry is suitable for studying en¬ 
tanglement and effects of interference which makes it an 
interesting device for quantum information processing. 
The flexibility of this setup makes it a convenient tool for 
investigating the transport properties of a TQD. Trans¬ 
port can be measured separately and simultaneously for 
the two double dot paths and be compared or combined 
to study the whole TQDs physics on the basis of the dou¬ 
ble dots. Also, transitions from double dot resonances in 
one path to configurations of all three dots in resonance 
can be studied in transport. In contrast to former pub¬ 
lished works [4] where one source and two drain leads 
were used, we now use one drain and two source leads. 
In this configuration of two-path transport the dot con¬ 
nected to the drain is shared by both paths (Figj^ (a)). 
The electrons from the different paths compete for the 
occupation of this dot. We analyze the role of interac¬ 
tions between the charge flowing through the two differ¬ 


ent paths by transport measurements. We observe, as a 
consequence of inter-channel Coulomb interaction, chan¬ 
nel blockade in transport. 
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FIG. 1. (a) Schematic of the TQD setup with capacitive and 
tunnel couplings, (b) AFM picture of the TQD sample with 
the in-plane gates Gi - G 4 and a QPC for charge measure¬ 
ments. The blue lines indicate the insulating barriers written 
by AFM. (c) Transport through path 1. Charging lines of 
dots A and B (solid lines) are observed. Charging of dot C 
(dotted line) is observed by a shift of the charging lines, (d) 
Transport through path 2. Charging lines of dot A and C 
(solid lines) can be seen and charging of dot B (dotted line) 
is observed by a shift of these charging lines. 

TQD Sample and Characterization. The measure¬ 
ments were performed on a lateral TQD made with 
local anodic oxidation by atomic force microscope 
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(AFM) on GaAs/AlGaAs-heterostmcture [2TH23]. A 
two-dimensional electron gas with an electron concen¬ 
tration of ne = 3,47 • is located in 33 nm 

depth below the surface. The dots A,B,G are arranged 
in a triangular geometry [4] with each dot placed next 
to the other two and one lead attached to each dot 
(Fig{^ (a)). Dots A and B and also A and C are tunnel 
coupled, dots B and C are capacitively coupled only. 
The source leads Si and S 2 are connected to dots B 
and C respectively and dot A is connected to the drain 
lead D. We have two transport paths: path 1 with dots 
A and B and path 2 with dots A and G. The sample 
has four in-plane gates G 1 -G 4 (Fig{^(b)) to control the 
potential of the dots, interdot and dot-lead couplings. 
A quantum point contact (QPC) sensitive to all three 
dots is placed next to dots B and C to perform charge 
measurements. The measurements were conducted 
in a dilution refrigerator. To measure the differential 
conductance of the two transport paths simultaneously 
but separately, a lock-in technique was used with ac 
voltages with two different frequencies fi = 83.3 Hz and 
f 2 = 18.3 Hz with Uac = 10 pV applied to Si and S 2 
respectively. The QPC was operated by applying a dc 
voltage to the source of the QPC, Sqpc, and measuring 
a dc current at the drain of the QPC, Dqpc. The QPC is 
tuned by the gate Gqpc- In our transport measurement 
range, the dots contain several ten electrons on the 
whole. The charging energies are Ech,A = 2 meV, Ech,B 
= 6 meV, Ech,c = 3 meV for dot A, B and C respectively. 

Charge Measurements. To characterize the device, the 
charging is studied by using the QPC as a detector. The 
derivative of the QPC current is plotted as a function 
of gate voltages Uqi and Uqs (Eigj^ with denoted 
charge configurations |Na,Nb,Nc), where Ni are the 
occupations of dots A,B,C. The electrons in the core of 
the dots are not included in Ni. The green lines indicate 
charging events, where one more electron is added to the 
system, pink lines indicate electron movement away from 
the detector. Charging lines with three different slopes, 
one slope for each dot, are visible in the measurement, as 
the slope depends on the capacitive coupling and thus on 
the distance between the dot and gates Gi and G 3 . The 
lines with the lowest slope belong to dot C as it is the 
least coupled to Uqi, the lines with intermediate slope 
to dot A and the lines with the largest slope to dot B as 
it is the least coupled to Uqs- Anticrossings of two dots 
in resonance are visible where two charging lines meet. 
At such a resonance two triple points (TPs) emerge 
where at each one three different charge configurations 
are degenerate. Charge reconfiguration lines connecting 
the TPs mark the charge transitions between the dots. 
Resonances between dot A and B (green circle), A and 
C (blue circle) and also between the only capacitively 
coupled dots B and C located in the two different paths 
(black circle) are observed [24]. When the double dot 
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FIG. 2. Charge measurement using the QPC. The green lines 
with different slopes point out the charging of the respective 
dot with one electron. Double dot resonances are marked by 
circles in green, black and blue, the triple dot resonance by a 
red circle. 

anticrossings coincide, the resonance condition for all 
three dots is fulfilled (red circle). We will focus on such 
a region in the transport measurements. 

Transport Measurements. To understand transport in 
this system, the differential conductance G is measured 
along path 1 and 2 simultaneously but separately, sweep¬ 
ing gate voltages \Jgi and Ugs (Figg (c),(d)). In doing 
so, the QPC is not in use. TPs with finite differential 
conductance are observed in both paths where the two 
dots are in resonance. In path 1 (Figg (c)) resonances 
of the dots A and B can be seen. Charging of dot C is 
observed by a shift of the charging lines of dot A and 
B, where dot C comes into resonance. Analogously, 
charging lines of dot A and C appear in path 2 (Figg 
(d)) and charging of dot B is detected indirectly by the 
shift of the charging lines of dot A and C. 

A triple dot resonance is formed where two double dot 
resonances coincide. In FigJ^(a) where we combine path 
1 and path 2 as observed in Figj^ (c),(d) we have three 
double dot resonances, A and B, A and C, B and C, in 
close vicinity to each other. One observes regions of high 
differential conductance in both paths but at different 
gate voltages Ugi and Ugs- Along the B charging line 
(path 1 ) we can identify two resonance lines where A is 
resonant with B. The whole triple dot physics can be 
observed. The different occupations of the states at the 
TP participating in transport are marked in Figj^ 

To understand the experimental results in more detail 
a theoretical model is developed which reproduces and 
explains the transport properties of the system. We fit 
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FIG. 3. (a) Combined color plot of differential conductance 
through path 1 and 2 with denoted charge conhgurations of 
the stability regions, (b) Calculated steady state occupation 
probabilities of each state i in the stability regions of 
the system (bottom). States 11,0,0) and |0,1,1) are a small 
region with low occupation probability. 


the transport simulations to the experimental data to 
extract the interdot tunnel couplings tab = 0.012 meV, 
^AC = 0.020 meV, the dot-lead tunnel couplings Fd 
= 0.008 meV, Fs^ = 0 . 003 , Fss = 0.006 meV and the 
electron temperature Tei = 300 mK. In the model we 
distinguish between particles coming from Si and S2 
in order to be as close as possible to the experimental 
conditions. 

Simulation. Using Master equation techniques (see 
Supplementary for more detail), we calculate the trans¬ 
port through the two paths and the stability regions of 
each state. 

In FigJ^ (b) we plot the numerical result of the steady 
state occupations pp. For small and large Uqi and Ugs 
the states |0,0,0) and are occupied respectively. 

Above |0, 0, 0), for larger values of Uqs, the state |0, 0,1) 
with one more electron in C becomes occupied and at 
the right hand side, for larger values of Uqi, the state 


| 0 , 1 , 0 ) with one more electron in B becomes occupied. 
Below |1,1,1), for smaller values of Uqs, the occupation 
of 11 ,1,0) with one less electron in C increases and at 
the left hand side, for smaller values of Uqi, the state 
| 1 , 0 , 1 ) with one less electron in B is occupied. All these 
regions obtained numerically correspond perfectly to 
the ones in FigJ^(a). The small regions of |1,0,0) and 
|0,1,1) connect the states with one electron and two 
electrons respectively. Each of these small regions, not 
really seen in FigJ^ (a), contain two TPs of path 1 and 
two TPs of path 2 . When two TPs coincide we have a 
quadruple point. 

With the information from the theoretical calculation of 
the occupations we determine the TPs present on the 
resonant lines in path 1 (FigJ^ (a)), (|0,0,0), |0,1,0), 
|1,0,0)), (|1,1,0), |1,0,0), |0,l,0))and(|0,0,l), |0,1,1), 
|1, 0,1)), (|1,1,1), |1, 0,1), |0,1,1)) where in the two last 
TP there is one more electron in C. In near vicinity of 
these TP there is high positive differential conductance 
in path 1 due to temperature broadening of the states. 
Thus they merge and form a vertical line of high differ¬ 
ential conductance. Similarly, for path 2 we have high 
positive differential conductance at the TPs (| 0 , 0 , 0 ), 
|0,0,1), |1,0,0)), (|1,0,1), |0,0,1), |1,0,0)) and, with 
one more electron in dot B, (| 0 , 1 , 0 ), | 0 , 1 , 1 ), | 1 , 1 , 0 )), 
(|1,1,1), |0,1,1), |1,1,0)). 

In Figj^ we plot the measured conductance of path 1 
and 2 separately ((a),(b)) as well as the results from 
the simulation ((c),(d)). In path 1 ( 2 ) we observe the 
splitting of the resonance between the dots A,B(A,C) 
due to the interaction with dot C(B). We also observe 
negative differential conductance in path 1 where path 
2 has high conductance (grey color in FigJ^ (a),(c)). In 
the following we will analyze and compare the transport 
features of the two paths in more detail. 

Channel Blockade. In FigJ^ we show a cut from 
Figi We observe that the resonance of path 2 splits 
into two (A,4) due to the interaction with the third 
dot present in path 1 . Path 1 gets in resonance at 
T and partially blocks the other path decreasing its 
conductance. This point is a quadruple point, four 
states of the two paths are coexisting in the same region 
of the stability diagram. The transport through path 2 
is stronger than through path 1 (tag > ^ab) thus when 
path 2 is in resonance (A,4) it totally blocks path 1 
decreasing its conductance even to negative values (more 
appreciable in Fig|^(a) and[^ (c))- 

This channel blockade is a consequence of Coulomb 
interaction between the charges flowing through the two 
transport channels which share dot A. When the bias 
voltage increases, the path with higher conductance 
increases its occupation in dot A blocking access to dot 
A from the other path and thus decreasing its transport 
and obtaining negative differential conductance. 

In Figj^ we identify the dominating and the blocked 
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FIG. 4. Differential conductance in experiment along path 1 
(a) and path 2 (b) and in simulation along path 1(c) and 
path 2 (d). The dotted line is the cut plotted in Figl^ 
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FIG. 5. Gut through the transport measurement and simula¬ 
tion of Figfor path 1 and 2 at Ug 3 = 103 mV. 


transport channels for each resonance with the 

information from the simulation. We show the initial 
and final states connected by two different transport 
paths, where in (A,4) one path blocks the other and in 
(t) both paths share the conductance. 

Blockade phenomena were previously studied for one 
dot attached to three leads (one drain and two sources) 
which contain some amount of up and down spins [ 25 ] . 
The spins of each path compete to occupy the dot 
blocking the access to the dot for the other spin by 
Coulomb interaction. Compared to this work for a single 
dot our TQD includes coherences between states of the 
two paths. Other papers [ 26 H 28 ] treat two transport 
paths just capacitively coupled where at some situations 
one of the paths blocks the transport through the other 
by Coulomb interaction. 

To observe the channel blockade in our experiment, it is 



FIG. 6. Transport mechanism for the peaks in Fig|^ In 
(A,^) path 2 blocks path 1 and in (t) both paths share the 
occupation of dot A. 


crucial to have two sources. If we switch the transport 
direction of both paths with the source lead connected to 
dot A and the two drain leads connected to dots B and 
C, the transport paths influence each other in a different 
way in the simulations. The electron flow splits at dot 
A in two paths with a probability that depends on the 
tunneling rates of path 1 (tab, Tb) and path 2 (tag, Tc). 

Conclusion. In summary, we have shown channel 
blockade in electronic transport through a TQD with 
two source leads. Coulomb interaction between electrons 
coming from the two sources gives rise to a blockade of 
transport through one path, when the other path has 
high conductance. The results show how interaction be¬ 
tween the joint transport paths is affecting the transport 
properties of the multi-terminal device and are a step 
towards a better understanding of transport properties 
in complex multi-dot systems. 

We are grateful to M. C. Rogge for producing the TQD 
sample. We acknowledge discussions with R. Sanchez 
and financial support from Spanish MICINN MAT 2014 - 
58241 -R 


kotzian@nano.uni-hannover.de 

[1] A. Vidan, R. Westervelt, M. Stopa, M. Hanson, and A. 
Gossard, Appl. Phys. Lett. 85, 3602 (2004) 

[2] L. Gaudreau, S. A. Studenikin, A. S. Sachrajda, P. Za- 
wadzki, A. Kam, J. Lapointe, M. Korkusinski, and P. 
Hawrylak, Phys. Rev. Lett. 97 , 036807 (2006). 

[3] D. Schroer, A. D. Greentree, L. Gaudreau, K. Eberl, L. G. 
L. Hollenberg, J. P. Kotthaus, and S. Ludwig, Phys. Rev. 
B 76 , 075306 (2007). 

[4] M. G. Rogge and R. J. Hang, Phys. Rev. B 77 , 193306 
(2008). 
































































5 


[5] B. Michaelis, C. Emary, and C. Beenakker, Europhys. 
Lett. 73 , 677 (2006) 

[6] C. Poltl, C. Emary, and T. Brandes, Phys. Rev. B 80 , 
115313 (2009). 

[71 M. Busl, R. Sanchez, and G. Platero, Phys. Rev. B 81 , 
121306 (2010). 

[8] C. Emary, Phys. Rev. B 76 , 245319 (2007). 

[9] R. Sanchez, E. Gallego-Marcos, and G. Platero, Phys. 
Rev. B 89, 161402(R) (2014). 

[10] E. Gallego-Marcos, R. Sanchez, and G. Platero, J. Appl. 
Phys. 117 , 112808 (2015). 

[11] M. Busl, G. Granger, L. Gaudreau, R. Sanchez, A. Kam, 
M. Pioro-Ladriere, S. A. Studenikin, P. Zawadzki, Z. R. 
Wasilewski, A. S. Sachrajda, and G. Platero, Nature Nan¬ 
otechnology 8, 261 (2013). 

[12] E. R. Braakman, P. Barthelemy, G. Reichl, W. Wegschei- 
der, and L. M. K. Vandersypen, Nature Nanotechnology 
8, 432 (2013). 

[13] R. Sanchez, G. Granger, L. Gaudreau, A. Kam, M. Pioro- 
Ladriere, S. A. Studenikin, P. Zawadzki, A. S. Sachrajda, 
and G. Platero, Phys. Rev. Lett. 112 , 176803 (2014). 

[14] L. Gaudreau, G. Granger, A. Kam, G. G. Aers, S. A. Stu¬ 
denikin, P. Zawadzki, M. Pioro-Ladriere, Z. R. Wasilewski, 
and A. S. Sachrajda, Nature Physics 8, 54 (2012). 

[15] D. DiVincenzo, D. Bacon, J. Kempe, G. Burkard, and K. 
Whaley, Nature 408 , 339 (2000). 


[16] M. Stopa, Phys. Rev. Lett. 88, 146802 (2002). 

[17] G. Lent, P. Tougaw, W. Porod, and G. Bernstein, Nan¬ 
otechnology 4 , 49 (1993). 

[18] M. G. Rogge and R. J. Hang, New Journal of Physics 11 , 
113037 (2009). 

[19] A. D. Greentree, J. H. Gole, A. R. Hamilton, and L. G. 
L. Hollenberg, Phys. Rev. B 70 , 235317 (2004). 

[20] H. Sasakura, S. Adachi, S. Muto, T. Usuki, and M. 
Takatsu, Semicond. Sci Technol. 19 , S409 (2004). 

[21] M. Ishii and K. Matsumoto, Jpn J. Appl. Phys. 34 , 1329 
(1995). 

[22] R. Held, T. Vancura, T. Heinzel, K. Ensslin, H. Holland, 
and W. Wegschneider, Appl. Phys. Lett. 73 , 2 (1998). 

[23] U. Keyser, H. Schumacher, U. Zeitler, R. Hang, and K. 
Eberl, Appl. Phys. Lett. 76 , 4 (2000). 

[24] J. M. Elzerman, R. Hanson, J. S. Greidanus, L. H. 
Willems van Beveren, S. De Eranceschi, L. M. K. Vander¬ 
sypen, S. Tarucha, and L. P. Kouwenhoven, Phys. Rev. B 
67 , 161308(R) (2003). 

[25] A. Gottet, W. Belzig, and G. Bruder, Phys. Rev. B 70 , 
115315 (2004). 

[26] R. Hussein and S. Kohler, Phys. Rev. B 89, 205424 
(2014). 

[27] G. Nietner, G. Schaller, G. Poltl, and T. Brandes, Phys. 
Rev. B 85 , 245431 (2012). 

[28] G. Schaller, G. KieBlich, and T. Brandes, Phys. Rev. B 
82 , 041303 (2010). 


Supplemental Material: Theoretical Model 


Here we discuss in detail the theory used to simulate the experiment. The total Hamiltonian of the system reads 
H = Hq + i^iead + ^int- The dot systcm (Ho) is a three site Anderson-like Hamiltonian [1]: + 

+ '^i<j where cj is the electron creation operator and fii the particle number operator of 

dot i. Ci {i = {A, B, C}) is the chemical potential of the dots, Tij the coherent interdot tunnel coupling and Vij 
the Coulomb interaction between the electrons in different dots. The reservoirs are modeled as a Fermi electron gas 
^lead = ^ikdlf^dik that has a constant temperature T and chemical potential {ii (/ = {Si, S2, D}). The interaction 
part of the Hamiltonian + H.c. couples the reservoirs and the dots with a hopping parameter 7/. 

The energy levels are tuned with the gate voltages present in the experiment (T/qi, Uq2^ Uq^)- 

The rates between the leads and the dots for incoming (+) and outgoing (—) electrons with respect to the dot 
system are given by Fermi’s golden rule = 27 r/h\^i\‘^f{/ai — e^) and = 27 r/h\ji\‘^[l — f{iii — e^)] where / 

is the Fermi distribution function. T/ = 27 r/fi| 7 /p is smaller than the interdot coupling thus we can apply the 
Born-Markov approximation [ 2 ] for the interaction of the system with the leads. From the Von Neumann equation 
= i/h[H^ which contains the full system time evolution, we trace over the baths degrees of freedom getting 
the reduce density matrix p{t) = Trieads^(^) [ 3 ] obtaining the master equation 

( 1 ) 

3 

pi it) is the occupation probability of the Tstate of the system and C is the Liouvillian superoperator that contains all 
the information about the system Hq and the jumping terms between the leads and the dots E^^/. As we just want to 
study the steady state properties of the system we solve the kernel of Eq. 0 to obtain the steady state occupations 
pSS _ Ker[£]. 

Taking the steady state occupations and the tunneling rates to/from the contacts we are able to calculate the current 



i,i=0 


( 2 ) 
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